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OPTIMAL FLOW THROUGH THE DISORDERED LATTICE^ 

By David Aldous 

University of California at Berkeley 

Consider routing traffic on the N x N torus, simultaneously be- 
tween all source-destination pairs, to minimize the cost ^ c(e)/^(e), 
where /(e) is the volume of flow across edge e and the c(e) form an 
i.i.d. random environment. We prove existence of a rescaled N —> oo 
limit constant for minimum cost, by comparison with an appropri- 
ate analogous problem about minimum-cost flows across a. M x M 
subsquare of the lattice. 

1. Introduction. In highly abstracted models of transportation or com- 
munication (e.g., roads, Internet) one is required (simultaneously for all 
source-destination pairs) to route a certain "volume of flow" from source to 
destination, and one seeks to minimize some notion of cost subject to some 
constraints (e.g., edge-capacities). In contrast to queueing theory, we shall 
regard flows as deterministic, but networks as random. A survey of such 
problems under various different models of random networks will be given 
elsewhere. In this paper we focus on the particular setting of disordered 7? 
and, after initial general discussion, study a specific model of costs. 

The mathematical structure of the lattice with an i.i.d. random envi- 
ronment (c(e) : e an edge of 1?) arises in many settings, such as first passage 
percolation [10, 14], disordered Ising models [4] and random walk in random 
environment [20]. Conceptually, the point of such disordered lattice models 
is to represent local spatial inhomogeneity. To study flow problems, consider 
a N X N square within Z^. To see how quantities scale suppose, for each pair 
{vi,V2) of vertices within the square, we send a flow of volume ajy from vi 
to V2 by spreading it evenly amongst minimal- length paths. Then the total 
volume of flow along a typical edge will be of order 

onxN'^x N/N'^ = N^aN, 
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where number of source-destination pairs = A^^; path length is of order 
N; number of edges is of order A'^^. We therefore scale by writing ajy = 
pN~^, and we interpret p as normalized traffic intensity. This normalization 
convention means that, in the optimal flow problems, we shall consider, 
flow volume along individual edges should be of order 1. To formulate one 
problem, suppose we impose capacity constraints: the total flow /(e) along 
each edge e cannot exceed cap(e), where the (cap(e) :e an edge of Z^) are 
i.i.d. with support bounded away from and infinity. It seems intuitively 
very plausible that there should be a critical value po such that 

P(feasible flow exists) ^1 as N ^ oo, p < po 

P(feasible flow exists) ^0 as N ^ oo, p > po- 

Moreover, it seems plausible that, if we also have i.i.d. costs (cost(e) :e an 
edge of Z^) representing cost per unit volume of flow across edge e, and if 
(when a feasible flow exists) we consider the cost C^v of a minimum-cost 
feasible flow, then we should have a limit function 

limiV~2ECiv = 5(p), P<Po- 

These are the kinds of results we would like to prove, but in fact (to avoid 
a specific technical issue discussed in Section 4.1) we replace "hard con- 
straints" cap(e) on the flow through an edge by "soft constraints" as follows. 
There is an i.i.d. environment (c(e) : e an edge of Z^) whose interpretation 
is now the cost of a total volume /(e) of flow across edge e equals c(e)/^(e). 
So the total cost of a flow is defined as 

(1) E^(^)/'(^) 

e 

summed over the edges e of the N x N square. Take normalized traffic in- 
tensity p = 1, and consider the cost Cat of a minimum-cost fiow (Cat is a 
random variable, the randomness arising only from the random environ- 
ment). Because flow volume across a typical edge should be of order 1, we 
expect a limit constant 

limiV-^ECAr = 7, 

depending only on the distribution of c(e). This is essentially what our main 
result (Theorem 1) establishes, with two further provisos. To exploit global 
stationarity, we work on the N x N torus instead of the square; and for 
technical reasons (Section 4.1) we impose a constant bound B on volume of 
flow across edges. 

Conceptually speaking, we are modeling congestion in the following sense. 
In a road network, the cost to users is the time taken to traverse a route. 
In the absence of congestion the cost per user does not depend on the flow 
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volume, so the total cost (all users combined) associated with an edge e is 
linear in /(e). If congestion causes speeds to decrease, it increases costs to 
all users of an edge, so the total cost is superlinear in /(e). The choice of 
/^(e) in (1) is just an arbitrary but mathematically convenient choice of 
superlinear function. 

Merely proving existence of a limit, as this paper does, is of course a 
rather modest achievement, so let us attempt several justifications: 

(i) Experience with other "disordered lattice" problems such as first 
passage percolation suggests it is very difficult to obtain explicit formulas 
for limit constants. 

(ii) A natural way to describe network behavior is via functions gNip) 
describing how some quantity varies with traffic intensity p. In the context of 
probability models for A'^-vertex networks, one would like to ensure existence 
of suitably renormalized N ^ oo limit functions g{p), and this paper is part 
of a project exploring methodology for such results. 

(iii) The standard technique for proving existence of limits in "spatial 
optimization" problems is subadditivity, which is readily applicable to, for 
example, traveling salesman type problems [18, 19], but which does not 
apply easily to our "flow" problems, because there is no simple way to relate 
the optimal flow on a 2N x 2N square to the optimal flows on the four 
N X N subsquares. Instead we use the local weak convergence methodology 
explained in Section 2.1. Our interest in this methodology arose from its 
use in optimization problems over locally tree-like random networks, where 
it becomes a reformulation of (special settings of) the cavity method [16] of 
statistical physics and allows explicit albeit nonrigorous calculation of limit 
constants. 

Plan of paper. Section 2 sets out notation, states the results carefully, 
and outlines the proof. The proofs themselves comprise Sections 3 and 4. 
We conclude in Section 5 with a wide-ranging discussion — about optimal 
flows in general, about local weak convergence methodology, about details 
of our particular model and about related work. 

2. Results and notation. Section 2.1 gives a rough verbal overview of 
methodology. We set up general notation in Section 2.2, enabling us to state 
our convergence result (Theorem 1) in Section 2.3. In Section 2.4 we set 
up notation specifically for dealing with flows across squares in the lattice, 
which permits us to state Theorem 2 encapsulating the methodological idea 
of relating flows within the N x N torus to flows across the M x M square. 
Section 2.5 then provides a more detailed outline of the proofs to follow. 
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2.1. Methodology. For M < N, it is obvious that an i.i.d. environment 
in a M X M subsquare of the N x N torus is distributed as an i.i.d. envi- 
ronment in a M X M subsquare of the lattice Z^. Our methodology rests 
on the less obvious idea that the — > oo limit of the global cost Cat on the 
torus (which involves route- lengths of order N), suitably normalized, can 
be identified with the M ^ oo limit of a certain quantity cm,b defined in 
terms of flows across the M x M square. We say across (from boundary 
to boundary) because as ^ oo the volume of flow originating at a ver- 
tex becomes negligible compared to the volume passing through the vertex. 
The argument has three parts. For a flow across a square, we call the joint 
distribution of entrance and exit points the transportation measure. 

Part 1. We consider — > oo limits of the empirical distribution, over 
all M X M subsquares, of the transportation measure for the optimal flow, 
jointly with the environment. 

Part 2. We use a concentration of measure argument to show that, given 
the empirical distribution of the transportation measure over a M x M 
square, for most realizations of the (c(e)) within the square the cost of 
optimal flow across the square is close to the mean cost. This allows us to 
take expectation over environments. Then convexity of the cost function 
(1) allows us to replace the empirical transportation measure by its mean 
measure. 

Part 3. The argument above leads to a definition (10) of cm,b which 
provides lower bounds on the N ^ oo limit optimal cost. To get the cor- 
responding upper bound, we need to construct flows on the N x N torus. 
Partition the torus into M x M squares. The definition of cm,b involves 
some particular transportation measure Qq across M x M squares. We con- 
struct flows by first constructing a "skeleton" as a Markov chain which steps 
from one boundary (between two squares) point to another boundary point; 
this skeleton does not depend on the realization of the environment. Within 
each square, use a minimum cost flow (which does depend on realization of 
environment) consistent with the transportation measure. 

2.2. General notation. Consider a graph G= {V,£), which for our pur- 
poses will usually be either the discrete N x N torus [0, — 1]^ (which we 
call T^) or the extended M x M square described in Section 2.4. A finite 
oriented path a = {vo,vi, . . . ,Vk) has a source or entrance ent{a) = vq and 
a destination or exit exi{a) = v^. For an edge e of G, write n{a,e) for the 
number of times that e occurs (in either direction) in a. Almost always we 
implicitly deal with self-avoiding paths, for which n{a,e) = or 1. A path- 
flow is a measure (by which we always mean a nonnegative finite measure) 
fj, on the space S of finite paths. The map 

o" — > (ent((T), exi((T)) 
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from S to V X V pushes /x forward to a measure on V x V which we call 
the transportation measure tra(//) by analogy with the classical mass trans- 
portation problem [17]. Call the edge- indexed collection 



the flow-volume flo(/x) associated with fj.. Thus, if tra(//) concentrated on a 
single element {v,v'), then flo(/i) would be a flow from source v to destination 
v' in the elementary sense of the max-flow min-cut theorem. But we use 
"flow" in the sense of multicommodity flow, and in informal discussion we 
do not distinguish carefully between path-flows and their associated flow- 
volumes. 

Consider a collection c = (c(e) : e an edge of £) of nonnegative real num- 
bers. Call c(e) a cost-factor and call c an environment. Given a flow-volume 
f and an environment c, define the cost of that flow in that environment to 
be 



To make our probability model, let the environment c be chosen i.i.d. from 
some distribution k with bounded support: for some < c* < oo, 



In what follows, E and P denote expectation and probability with respect 
to the random environment. 

Call a measure on a finite set B constant if 6{b) is constant for all 6 € .B; 
we reserve uniform for the case of a probability measure. 

Write M\N for 'W is a multiple of M." 

2.3. The main theorem. Now take the graph to be the discrete torus 
Tjy = [0, — 1]^, and where helpful append a subscript -(^v) to the notation 
above, so that, for instance, 5^(7v) is the set of finite paths in T^. As explained 
informally in the Introduction, we want to consider flows with volume N~'^ 
between every source-destination pair. In the terminology above, define a 
standardized global flow f to be the flow-volume of some path-flow ^ for 
which tra(/i) is the constant measure 




(2) 




(3) 






6 



D. ALDOUS 



For a fixed environment c, consider the minimum cost subject also to this 
capacity constraint: 

cost(jv)(c, -B) := inf{cost(f , c) :f a standardized global flow satisfying (4)}. 
Now make the environment random as at (3). 

Theorem 1. There exists a constant ^{k,B) such that 
lim iV~^Ecost(Ar)(c,S) = 7(k;,S). 

The corresponding SLLN then holds by a concentration inequality; see 
Section 5. 

2.4. Flows across subsquares: notation. We now develop notation for use 
in the proofs. 

Figure 1 (middle right) shows the M x M square [0,M-1]2 (regarded as 
a subgraph of the lattice Z^) together with half of each edge from this square 
to its complement. Create artificial boundary points b at ends of these half- 
edges; each such b is of the form (— ^, j) or (M — |, j) or {i, — i) or (i, M — i). 
Write for the set of all internal vertices. Write Boum for the set of 
4M boundary points, and write Bm = B%j U Boua/. Write £m for the edge- 
set (internal edges and boundary half-edges). Call {Bm,£m) the extended 
M X M square. We use the general notation of Section 2.2 for the graph 
{Bm,Sm), writing a subscript -m where helpful. So T,m denotes the set of 
finite paths in this graph, and a path-flow /i has a transportation measure 
tra(^) on Bm x Bm and a flow-volume f = flo(/i) on the space Tm C [0, oo)^*^ 
of possible flow- volumes. Write c = (c(e), e € £m) for an environment of cost- 
factors on the edges of the extended M x M square; write Cm for the set 
of all such environments c. For a flow- volume f and an environment c, we 
write 

(5) costM(f,c)= c(e)/'(e). 

The preceding just copied general notation; we now come to new defini- 
tions specific to the extended M x M square. Partition the boundary points 
BouAf into four subsets 

BouM = Boua7 U Bou^f ^""^ U Bou^ff U Bou^^f * 

in the way implied by the notation (Figure 1, bottom right). For b G Boumj 
let b^^^'^'^^ G BouAf be the boundary point obtained by reflecting b top-to- 
bottom or left-to-right; in particular, ft^'^fl'^ct _ ^ j^q(^ (M, M) (Figure 1, bot- 
tom right). 
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N X N torus 



vertices 



extended M x M square 



vertices 
Bm 



vertices 



Boy;/. 



vertices 
Boirif 



Fig. 1. Top left: the N x N torus (here N = 20), in which the edges wrap from left side 
to right side and from top to bottom. Middle left: the natural partition of the N x N torus 
into M X M subsquares. Here M = 5. Bottom left: the boundary points Boujv.a/ for the 
partition. Top right: the "ordinary" M x M square [0,M - if. Middle right: the M x M 
square with its boundary points. This is what we call "the extended M x M square" in the 
text. Bottom right: The set Bouj\/ of boundary points of the extended M x M square, and 
the mapb^^"^""'-. 
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Given a transportation measure Q on Bou m x Bou m , write Qcnt and Qcxi 
for its marginals on Boum (later we use the same notation for transportation 
measures on other grids). Define Qm to be the set of transportation measures 
Q such that 

(6) the push-forward of Qexi under b —> b^^^^'^^ equals Qent- 

For Q € Qm , define 



1 



(7) drift(Q) = ^EE(^2 - h)Q{bi,h 



So drift((5) is a point in M^. Taking this point modulo (1,1) gives a point 
in the continuous torus = [0, 1)^: 

drifti(Q) := drift(Q) mod(l, 1). 

Finally, define Qo to be isotropic if it has a representation as a mixture 



(8) Qo=/ Q^{dQ), 

JQm 

where ^ is a probability distribution on Qm whose push- forward under the 
map Q drifti((5) is the uniform probability distribution on [0, 1)^. Write 
Qf^ for the set of isotropic transportation measures. Because property (6) 
is preserved under mixtures, Qff C Qm- 

Given a transportation measure Q on Boua/ x Bou a/, and an environ- 
ment c, define 



(9) costjvf bIQ, c) = inf < costj./ (flo(/.i), c) : tiahi) = Q, maxflo(/u)(e) < B >. 

^l [ e ) 

We now arrive at the central definition: 

(10) CM,B :=inf{EcostM,B(Q,c) :Q e Q'H}. 

Here E is expectation w.r.t. the random environment, with the following 
convention. Cost-factors c(e) are i.i.d. (k) as e runs over internal edges and 
over edges to boundary points on the bottom and left sides; c(e) = on other 
edges to boundary points. Obviously, this convention is designed so that 
when the N x N torus is partitioned into M x M squares (Figure 1, middle 
left), each edge of the torus is assigned to a unique square. In bounding costs 
we may take the number of edges \£m\ as 2M^. 
Here is the promised elaboration of Theorem 1. 

Theorem 2. For B > 1/4, there is a limit constant 



And 



7(k,S)= lim M-^CM,B- 



lim iV"^Ecost(7v)(c,S) =7(k,S). 
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2.5. Outline of proof . Here we expand the rough description of method- 
ology (Section 2.1) into a more detailed outline of the proof. This outline 
will ignore the constraint /(e) < B, which enters only at a technical level, 
sometimes as a convenience and sometimes as an inconvenience. 

Getting an upper bound on the limit lim7v-+oo N^'^K cost(7v) (c, B) requires 
a construction of a standardized global flow. The definition (10) of cm,b 
involves some particular transportation measure Q € attaining the in- 
fimum. This Q is a positive matrix on the boundary vertices of a standard 
extended M x M square, but can be used in a natural way to define a 
Markov chain on the boundary points Bou7v,a/ of the partition of the N x N 
torus into M x M squares: see Figure 2. Condition (6) is the essential "com- 
patibility" condition making this construction work. Running the chain for 
a suitable number of steps with a suitable starting distribution on Boutv.a/ 
gives a joint distribution of starting and ending points which is approx- 
imately uniform on Bouat^a-/ x Bou]y,M- The definition (8) of isotropic is 
what gives this "approximately uniform" property (Proposition 4). Making 
the walk into a standardized global flow requires appending length o{N) 
path-segments at the ends to make the transportation measure become ex- 
actly uniform, and expanding a step between two points on the boundary 
of an M X M square into a path across the square. But the construction 
makes the transportation measure across a square be the original Q, so to 
"expand a step" we simply use the routing attaining the minimum in the 
definition (9) of costj\/,B(Q, c); see Proposition 5. This completes the con- 
struction of a standardized global flow. The cost associated with the end 
o(A^)-length segments is negligible, and the cost associated with flow across 
a typical M x M square in the torus is cm,b, leading to the desired upper 
bound (Proposition 17). 

The lower bound uses more abstract methods. Consider the optimal flow 
within the N x N torus, then consider this flow within a randomly-positioned 
M X M square, and then consider subsequential N ^ oo weak limits for fixed 
M. What can we say about such a limit distribution? At first sight we cannot 
say anything explicit, because we do not know anything explicit about the 
optimal flow on the N x N torus. But we can collect some properties, as 
follows: 

(i) The environment c = (c(e)) within the M x M square follows the 
original i.i.d. model (3). 

(ii) Because the volume of flow within the torus originating or ending at 
a given vertex is N~'^, in the present N ^ oo limit all the flow is across the 
M X M square from boundary to boundary. 

(iii) There is some random transportation measure Q giving the joint 
distribution of entrance and exit points of flow. Note Q will be dependent 
on c. 
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Fig. 2. Left side: 4 steps (5^0, ^i, ^2, Va, Ki) of a walk on the skeleton graph. Right side: 
the projection walk (projy ^,j(yo),prOjv,M(^i).PrOjv,A/(^2),prOjv_j,.j(y3),prOjv,M(^4)) on the 
boundary of the extended M x M square. The state of the projection walk is the relative 
position of entry to the next square for the skeleton walk. 

(iv) Because a global optimum is a local optimum, given the environment 
c and the transportation measure Q, the flow- volumes f = (/(e)) within the 
M X M square minimize the local cost 

costM(f,c)= ^ c(e)/^(e) 

subject to the given transportation measure. 

(v) The fact that the optimal flow on the torus was some standardized 
global flow constrains the induced flows across a typical M x M square, and 
what it implies is (roughly, and in part) that EQ G Qj^f • 

(vi) The mean cost on the torus relates to the mean cost associated with 
our typical M x M square. 

Lemma 20 formalizes these assertions in the finite context. Letting 
— > oo gives a lower bound for limjv^oo A''~^Ecost(Ar)(c, B) involving 

(11) EinfjcostAf (f,c); f has transportation measure Q}. 

Lemma 21 states this bound precisely. Of course, we do not know explicitly 
what is the distribution of Q, so we need to lower bound the quantity (11) 
by appealing to constraints Q must satisfy. To do so, the key ingredients are 
the two very easy facts referenced below. Consider a nonrandom Q. Then by 
an easy concentration inequality (28), the random variable costM, BiQ,c) is 
close to its expected value E cost m,b(Q5 c). So if Q were random but took val- 
ues in a small set (Qj) of possible values, then the expectation in (11) would 
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be approximately the appropriate weighted average of EcostA,/,_B(Qj, c). Be- 
cause (Lemma 18) Q — >■ EcostM,s(Q) c) is convex, we could then lower bound 
the quantity (11) by ¥.cosiM,B{Q-,c) for Q = EQ € Qf°, which by definition 
is lower bounded by cm,b, giving the desired lower bound in Theorem 2. 

So the remaining issue is to show that the quantity (11) is not much 
changed by replacing Q by a "quantized" version taking values in some 
large finite set of "smoothed" transportation measures, which we define as 
follows. Set M = KL and partition the boundary of the M x M square into 
blocks of length K. Require a smoothed transportation measure Q' to have 
the property that Q'(6i,62) depends only on the blocks containing bi and 
62, and that Q'{bi,b2) be an integer multiple of some small constant. Then 
it is enough to show that, in an arbitrary environment c, the cost of the 
cheapest fiow f with transportation measure Q is not much less than the 
cost of the cheapest fiow f with smoothed transportation measure Q' . This 
scheme turns out to be technically complicated to implement (see Section 
4.3), but a variant is carried though in Lemma 19. 

3. Constructing global flows from local flows. In this section we show 
how to use flows across the extended M x M square to construct stan- 
dardized global fiows on the N x N torus. This leads to the upper bound 
of Proposition 17. The first stage of the argument leads to a "clean" con- 
struction of a global fiow (Proposition 5) for which the source-destination 
distribution is approximately constant. We then need to "patch" to make 
the source-destination distribution be exactly constant; this is done using 
some of the elementary constructions collected in Section 3.2. 

3.1. The basic construction. Fix M\N. Partition the N x N torus into 
M X M squares in the natural way (Figure 1, middle left). Insert bound- 
ary points midway along edges linking different squares. This creates a set 
Bouat^m of 2N'^ /M boundary points. Define the skeleton graph (Boun,m,£n,m) 
to have an edge between each pair of boundary points which are boundary 
points of some common square S. Consider a measure z/ on paths (xj) in 
the skeleton graph. For any square S with boundary Bou(5), and for each 
(61,62) GBoU(5) X BoU(5), write 

(12) u^s){bi,b2) =J2'^{i^o, xi,X2, ...)■■ {xi, Xi+i) = (61,62)} 

i 

for the mean number of times the path steps from 61 to 62. We often discuss 
a typical M x M square as if it were the standard M x M square by saying 
"up to translation." 

We define some properties for a measure 9 on Bouat^m x Bouat^m with 
marginals Bent and ^cxi- 
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Properties 3. (i) ^ has total mass N. 

(i) 6 is invariant under translation of the torus by (M, 0) or (0,M). 

(ii) ^ent = ^oxi ■ 

Recall that Qm is a set of transportation measures on the boundary of 
the extended M x M square. 

Proposition 4. Given Q € Qm satisfying an irreducihility property 
(14), there is a measure 9n,m on Bouat^a/ x Bouat^m satisfying Properties 3 
and the following two properties: 

(i) There is a measure v on paths in the skeleton graph whose transporta- 
tion measure tra(z^) = 9n,m o,f^d such that, for each square S , the measure 
1/(5) at (12) equals Q (up to translation). 

(ii) Let 9n,m be the push-forward of 9n,m under the map {z^,z'^) —>■ 
{j^z^j-^z"^) from Bouat^a/ x Bou7v,m to x (where is the continu- 
ous torus [0, 1)^J. Then as N ^ oo, there is weak convergence of N~^9j\f^M 
to the distribution of {U,U + drifti((5)), where U has uniform probability 
distribution on . 

Proof. Use Q to define a transition matrix Q on Boua/ via 

Q(6,6') = Q(^fe')/Qcnt(6). 

This is a (Markov) transition matrix because Qent{b) = Z^b' Q{b,b'). The set 
Bouat^a,/ of inter-square boundary points b can be expanded to a set Bou^ ^ 
with elements 

(6,^) and (6,^); or (6, j) and (6, J,) 

indicating a direction along an inter-square edge. Use the transition matrix 
Q to define a transition matrix Q^{-, •) on Bou^ ^ as follows. For a state, say, 
(60,^)5 the arrow points into some square 5o = [xo,a;o + M — 1] x [yo,yo + 
M — 1], and 60 is of the form bo = (xo,yo) + &o some 6q G Boua/. Any 
bl E Boua/ determines a point 61 = (xo,yo) + ^1 € BouAr,A/- Define 

gi((6o,-),(6i,tt)) = Q(6S,6D, 

where [| is the direction of arrow pointing out of Sq at 61. 

Using the defining property (6) for Q G Qm, one can check that a station- 
ary distribution vr on Bou^ for is 



(13) 



T^{{xo,yo) + bl4) =Qcnt{b*i)/q, 
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where as above fj is the direction of arrow pointing out of Sq at 61, and where 
the normahzation constant is g' = (^)^Q(Boujvf x Boum). Define 

q__ NQ(BouM X BouAf) 

and suppose first that t is an integer. Define a probabihty measure iy on 
paths in the skeleton graph by: 

• picking an initial state in Bou^^^j from distribution ir, 

• running the Q^-chain for t steps, 

• deleting arrow-labels. 

Consider a typical element {bo, —^) of Bou^ jyj with bo = {xq, uq) -|-6q for some 
6q G Bou7\/. Recalling definition (12) and using stationarity in the first line 
below, 

P(5)(6o5 ^1) = ti>{ paths with first step 60 &i} 
= t7r(6o,^)Q'((6o,^),(&i,W) 

_^ Qent{hl)Q{bl,bl) 
q Qcnt(bS) 

and the same identity holds for general (60 ,bi). Now define the normalized 
measure i^(-) ■.= Ni'{-). Then 

fN 

by definition of t. Define 6n,m = tra(i^). Properties 3 follow from the con- 
stancy properties of the stationary distribution vr at (13), and assertion (i) 
of Proposition 4 is immediate from the construction. 

Now note that if t is not an integer, then we can apply this construction 
to [t\ and to [t] , and mix over these two cases. 

It remains to prove assertion (ii). Note that the stationary vr at (13), 
projected down to Bouat^m, is asymptotically uniform as ^ 00, in the 
sense that its push- forward under the map z — > N~^z converges weakly to 
the uniform distribution on the continuous torus T^. 

Note that Q defines a Markov transition matrix Q* on Boum via 

Q*(6o,6i) = g(6o,6r*^^^*), 

where fei ^reflect jg ^j^g "refiect top boundary with bottom boundary, refiect 
left boundary with right boundary" procedure above (6). 
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Define a projection map pro^Y j^j : Bou^ — > Boum as "take modulo (M, M), 
choosing the boundary point so that the arrow points into the standard 
M X M square." So, for instance, 

pro^,j,^((3Af + i,2M- i), j) = (i, -i) 

pro^,M((3M + i,2M- i), [) = (i, M - i). 

Here is the key idea in the construction, illustrated in Figure 2. Let (Kj, s = 
0,1,2,...) be the Q^-chain, or, more precisely, the chain defined in the same 
way over the lattice Z^, so that its values modulo {N,N) are the Q^-chain. 
Then it is straightforward to verify 

(Xs ■■= pro^ s = 0, 1, 2, . . .) is the (5*-chain. 

So the displacement It — of the Q^-chain over T steps can be written as 

T 

Yt-Yo = Y. g{Xs-i,Xs); g{bo,h) = bf''^'' - bo- 

s=l 

Note the right-hand side does not involve A^. Suppose first 
(14) the Q*-cham is irreducible. 

Then by the strong law of large numbers for additive functionals of a finite 
Markov chain, 

T 

s=l b,b'GBouM 

where vr* is the stationary distribution of Q*. In terms of t = 
steps of the Q^-chain (Yg), this says that as — > oo, 

g(3(BouA/ X BOUA/) . 
N \Yt- Yo) m probability. 

To calculate g, note that 7r*(-) is proportional to (5cnt(')- So 
EbM'{b'-b) Q{b,b') 



Q(Bouj\/ X BouA/) 
M2 drift(Q) 



Q{BouM X BouA/) 

and so 

(15) iV"^ {Yt - Yo) drift(Q) in probability. 
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Looking at assertion (ii) of Proposition 4, we noted earlier that the first 
marginal of N^^9N^]\f converges weakly to dist(C/), and (15) now shows 
that the conditional distribution converges to the unit mass at drifti((5), 
completing the proof. □ 

Recall the definition (10) of cm^b- By compactness, the inf is attained, 
so we have 

CM,B = IE costM,s (Qo , c) , 

where 

(16) Qo=/ Q^{dQ) 

Jqm 

for a certain probability measure "0 on Qm- 

Proposition 5. Given an environment c on the N x N torus, we can 
define path-flows fiN, Aii'lc) with associated flow-volumes f(-|c) on the torus 
such that: 

(i) maxg /(e|c) < B. 

(ii) Ecost(yv)(f(-|c),c) = {^fcM,B- 

(iii) tra(/.fAr^A/(-|c)) is a measure pN,M on Bouat^a/ x Bouat^m which does 
not depend on c, and which satisfies Properties 3. 

(iv) Let pN,M be the push-forward of pN,M under the map {z^,z^) — > 
{^z^^j^z^). Then as N ^ oo with M fixed, N~^p]\f^M converges weakly to 
the uniform probability distribution on x T^. 

Proof. Deferring the issue of irreducibility in its hypothesis, Proposi- 
tion 4 associates with each Q £ Qm a measure on skeleton-paths. Define 
the mixture 

Jqm 

corresponding to (16). Then pn,m '■= tra(i^) satisfies Properties 3 by Propo- 
sition 4 and the fact that these properties are preserved under mixtures. 
Result (iv) follows from Proposition 4(ii) and the definition of isotropic [un- 
der tjj, the distribution of drifti((5) is uniform on the continuous torus T^]. 
The result that 

(17) iy(s) = Qo 

follows from Proposition 4(i) and the mixture construction. 
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Now consider the extended Af x Af square. By definition of cm,b, there 
exist path-flows //"(-Ic) across the square (depending on the environment c), 
with flow- volumes f'^(-|c), such that 

(18) max/°(e|c)<B Vc, 

(19) tra(/i°(-|c)) = Qo Vc, 

(20) EcostM(f°(-|c),c)=CA/,B. 

For b,b' G Boum, write /i|^{j,(-|c) for the flow ;u'^(-|c) restricted to paths with 
entrance b and exit b' , and normahzed to have mass 1. 

Recalhng the construction of in Proposition 4, we see that u can be 
constructed in 5 steps: 

• pick Q from distribution ■0, 

• pick an initial state xq in Bou^^j- from the stationary distribution of Q, 

• run the Q^-chain for t steps, 

• delete arrow-labels to get a path (xo,xi, ■ ■ ■ ,xt) in the skeleton graph, 

• take the probability distribution of the skeleton-paths thus defined, and 
multiply the measure by N. 

This gives measure ly. But we can augment the construction in the natural 
way. Let c be an environment on the N x N torus. In each step of the skeleton 
walk, the successive points {xi-i,Xi) are in the boundary of some subsquare 
Si. Write c*^*) for the restriction of c to Si. Expand each step (xj_i, Xj) into a 
path {xi-i = yQ,yi, . . . ,ym = Xi) chosen from distribution M2i_i,a;,('l'^''*^) (up 
to translation) . Augmenting the construction in this manner gives a measure 
/^Ar,A'/('|c) on paths in the N x N torus. It is easy to check [the key point 
being that the same transportation measure Qq appears in (17) and in (19)] 
that the restriction of these path-fiows to a subsquare Si gives path-fiows 
on Si agreeing (up to translation) with the path-flows /i'^(-|c) across the 
standard M x M square. So assertions (i) and (ii) follow from (18) and (20). 

Finally we outline the technical issue of handling the irreducibility con- 
dition in Proposition 4. In the mixture representation (16) let T be the 
push-forward of -0 under the map Q — drift(Q). It is not hard to show 
[cf. Lemma 20(vi)] that J^2 \\u\\iT(du) < oo. Make some particular choice 
of irreducible transportation measures Q^"^ G Qa/ with drift(Q(")) = u and 
construct the mixture 

Q*= [ Q(d"ft{Q)) ^(dQ). 

Then eQq will be a feasible transportation measure, for sufficiently small e. 
So we can repeat the proof above with the mixture (16) replaced by 

(1 - e)Qo + eQl= j ((1 - e)Q + eQ^^^^ft^Q))) V'((iQ) 
-IQm 
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because now each (1 — £)Q + eQ(<^^'^*('3)) is irreducible. This gives the same 
conclusion except that in (ii) cm,b is replaced by cm,b{£) := IEcostM,B((l — 
e)Qo + ^Qo)C). Because cm,b{£) — > cm,b as e | 0, this conclusion suffices for 
our later needs. □ 

3.2. Patching and smoothing flows. Here we collect an assortment of 
simple lemmas, mostly involving flows whose construction does not depend 
on the environment. Some are used in the next section in proving the upper 
bound, the others in Section 4.6 in proving the lower bound. 

Write f""^' for the flow-volume on the N x N torus associated with the 
standardized global flow which puts equal weight on all minimum-length 
paths between entrance and exit. 

Lemma 6. /''"(e) = [N^ /2\/{2N'^) < \ for each edge e. 

Proof. Write ||f — for the graph-distance (i.e., shortest path length) 
between vertices v and w on the torus. By symmetry, /""'(e) does not de- 
pend on e, so 

2iVV"'"(e)=E/"'(e) 

e 

V w 

= Ar"i^||^_(0,0)||i 

V 

N-l 

= 2 ^ mm{i, N — i) 

i=l 

= [iVV2j • □ 

The next lemma has a similar proof, which we omit. 

Lemma 7. Let U be a uniform random vertex in the N x N torus, and 
let U(^i) be a uniform random point in the L x L square centered at U , where 
1 < L < N is odd. Then there exists a path-flow fi such that 

tra(/i) is the distribution of {U,U(^x,)) 



maxflo(^)(e) < 
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Lemma 8. Let p be a measure on the boundary Bou^/ of the extended 
M X M square, and let p be the uniform probability distribution on the in- 
ternal vertices Bou^^. There exists a measure p on path-space such that 

tra(/i) = p X p, 
maxflo(/u)(e) < 2maxp(5). 



Proof. Use the path from a boundary vertex b to an internal vertex a 
which makes at most one turn. Check that, for any e, the number of pairs 
(6, a) for which the path traverses e is at most 2M^. So 

/(e)<2M2.maxp(6).^. □ 

Lemma 9. Let Q be a measure on Boujv/ x Boua/. There exists a measure 
p on path-space Sm such that 

tra(/i) = Q, 

maxflo(/_f)(e) < 2( maxQcxi(^) + max(5cnt(^) I • 

Proof. For each pair (6,6') € Boum x Boum, route flow via a uniform 
internal vertex, and use Lemma 8. □ 



Lemma 10. Let 9 be a measure on the bottom boundary points Bou^**"™ 
of the extended K x K square, and let 9 be the uniform probability measure 
on the top boundary points Bou^'^. Then there exists a path- flow p such 
that: 

(i) iTa.{p) = 9x9; 

(ii) maxe flo(//)(e) < maxft 0(6). 



Remark. We use K instead of M here for consistency with its applica- 
tion later. 



Proof of Lemma 10. By scaling, we may suppose maxb6'(6) = 1. Then 
it is enough to consider the special case 9{b) = 1 V6 G Bou^""""^ and delete 
surplus flow. In this special case, we route flow from (x, — |) G Bou^*^*"™ to 
(y, if — ^) € Bou^^ as follows. If x = y, we route straight upward. Otherwise 
set d=\y — x\ — 1 and route from (x, —\), up to (x, d), across to {y, d), up 
to iy,K-^). 



FLOW THROUGH DISORDERED LATTICE 19 

It is clear that the flow on each vertical edge equals 1. For the horizontal 
edge {i — l,d) — {i,d), the flow equals 

2 

— \{{x, y):0<x<i-l, i<y<K-l,y = x + d + l}\ 
K 

2 

< — max(i, ii' — i) < 1. □ 
K 

The next three lemmas compare the costs of different flows. Note that the 
same device [the factorization in (21)] is used each time. 

Lemma 11. For\<Bi<B2, 

Ecost(7v)(c, S2) < Ecost(7v)(c, 

* ,r2-^2(-B2 - ^l) 



<EcOSt(7V)(c,S2) +C*A^^ 



B2 - 1/4 



Proof. The left inequality is immediate. Fix c with maxec(e) < c*. 
Let f^^^ be a flow volume for which max^ /(2)(e) < B2 and cost(f(2),c) = 
cost(7v)(c, B2) and let P"^' be the uniform flow from Lemma 6. Deflne f^^^ = 
^funi ^ (1 _ ^)f{2)^ ^ ^ jg ti^g solution of 

Bi = \X + B2{l-\). 
Then f is a standardized global flow with maXg /(i)(e) <5i, and 

cost(fW,c) -cost(f(2),c) 

(21) 

<c*{2B2)Y.{\\) 

e 

= c*B2N^\ 
establishing the lemma. □ 

The next lemma relates almost-isotropic transportation measures to isotropic 
ones. Let be the path-flow across the extended M x M square which 
assigns weight 1 to each left-to-right horizontal path; let be its trans- 
portation measure. Deflne similarly n^. 



Lemma 12. Suppose Q G Qm and 
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Write e = ei + £2 + £3 + £4 . Then for B' > B + e, 

EcostM,ij(Q,c) > cm,b' - 4:C*B'eM'^. 

Proof. Given c, let /x(-|c) be a path-flow attaining the inf in the defi- 
nition (9) of costjvf,B(Q,c). Consider 

whose transportation measure is 

Q':=Q + eiQ^ + £2(3^ + esQ^ + e^Q^ G Qj^°. 

Write f and f for flo(^(-|c)) and flo(;u'(-|c)). Note maxe/(e) < B + e < B' 
and 

costM(f',c) - costM(f,c) = ^c(e)(/'(e) + /(e))(/'(e) - /(e)) 

e 

<maxc(e) x |£^7v/|(2i?')£. 

So 

CM,B' < ]EcOStM,B(Q',c) 

<EcostA/(flo(^'(-|c)),c) 
< EcostM(flo(^(-|c)),c) + 2c*B'\SM\e 
= KcostM,BiQ,c) + 4:C*B'M'^e, 
where we have used \£m\ = 2M2. □ 

The next lemma relates the cost of an almost-standardized global flow to 
the cost of a standardized global flow. 

Lemma 13. Let f = fi -|- £2 be the flow-volume in a standardized global 
flow on the N x N torus. Let < 6 < Bq — ^. Suppose 

max/i(e) < Bq, 
max /2(e) < 5. 

Then there exists a standardized global flow whose flow-volume f satisfies 
max /(e) < Bq — 5, 

rn . . 85N^B^ 

cost(f, c) < cost(fi, c) + — — — maxc(e) 

-DO — 1/4 e 

for every environment c. 
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Proof. Let f"°' be the "minimum-length" flow-volume from Lemma 6, 
so maxg /""(e) < 1/4. For < A < 1, consider the mixture ix = {I — A)(fi -|- 
£2) + Af*^"' arising from the corresponding mixture of path-flows. Choosing 
A ~ s+Bo-i/A which is the solution of 

Bo-d={l-\){Bo + 6) + lX 

ensures that maxg fx{e) < Bq — 6. For any e, 

flie) - f!{e) = {fx{e) + h{e)){fx{e) - /i(e)) 
<max(0,2So(/A(e)-/i(e))) 
<2Bo{S + Xr\e)) 

- Bo -1/4' 

bumming over the 2iV2 ed ges e gives the stated bound for cost(f,c). □ 

Lemma 14. Identify the edges of the N x N torus with a subset of the 
edges of the (N -|- 1) x (N + 1) torus in the natural way, identifying edge 
(0,i) — (A^ — l,i) with edge (0,i) — {N,i). Let c' be an environment on the 
(N -|- 1) X (N + 1) torus with maxec'(e) < c* , and let c be the environment 
on the N X N torus induced by the identification. Let f be a standardized 
global flow on the N x N torus such that maxe/(e) < B. Then there ex- 
ists a standardized global flow f on the {N -|- 1) x {N + 1) torus such that 
maxe /'(e) < B and 

cost(f' ,c') <cost(f,c) +{N +l)G{c*,B), 

where G{c*,B) depends only on c* and B. 

Proof. We construct a flow in three pieces. A path in the N x N torus 
induces a path in the (iV -|- 1) x (N + 1) torus in a natural way, where 
necessarily identifying an edge such as (0,i) — {N — l,i) in the former with 
the two edges {0,i) — {N,i) — (N — l,i) in the latter. So the flow f induces 
a flow fi with associated path-flow /ii on the {N -|- 1) x (A^ -|- 1) torus such 
that, after renormalizing to ensure 

tTa{fii){vi,V2) = , _^ .3 , vi,V2 e[0,N - 1]2, 
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we have 

(22) max/i(e)<i?(^^^j , 

(23) E <^)fhe) = (-^)'cost(f,c), 

where ^'at denotes the subset of the edges of the {N + 1) x (iV + 1) torus 
identified with the edges of the N x N torus. 
For the second piece, consider the vertex set 

A := {{i, N):0<i<N-l}U {{N, i) -.0 <i < N - 1}. 

We shall define a path-flow fi2 such that 

(24) tra(^2)(^^l,^^2) = (iV + l)"^ {vi,V2) e A x [0, N - if . 

We can do this using Lemma 8 with M = N, because the boundary points 
of the extended N x N square can be mapped via a two-to-one map to A. 
Lemma 8 with p uniform, under this mapping, gives a path-flow /i2 such 
that 



tiSi{fL2)ivuV2) = 2x ^ X {vi,V2)eAx [0, iV - 1]^ 

(25) maxflo(/i2)(e)<4x^. 

Renormalization gives /X2 satisfying (24) and 

2N^ 



maxflo(^2)(e) < 



{N + lf 

For the third piece we want a path-How fi-^ such that 

tra(M3)(^i,^2)= ^^l^p , 

{vi,V2)eAx AuAx {{N,N)}U{{N,N)} x A. 

We can construct such a flow using only the edges of £(n+i) \£{n) which 
are parallel to the boundary of the N x N square, that is, of the form 
(A^, i) — {N, i -|- 1) and (z, N) — (i + l, N). It is easy to construct such a flow 
with 

maxflo(/i3)(e) < ^^qjy- 
These edges were not used by the previous flows. 
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Now superimpose the path-flow ^2 with its reversal (whose transportation 
measure is supported on [0, — 1]^ x A) and with ;U3, and write £2 for the 
resulting flow-volume. So we have constructed flow-volumes fi and £2 on the 
(iV -I- 1) X (A^ 1) torus such that 

fi -|- £2 is a standardized global flow 

max/i(e)<S by (22), 

max/2(e)<^^ by (25), 

cost(fi,c') < cost(f,c) 2iVc*S^ 

the latter by (23) and the crude bound for the doubled edges. Applying 
Lemma 13 to A^ -|- 1 and with 5 = establishes the result. □ 

Finally we state a result, somewhat analogous to the previous lemma, 
relating cm+i,b to cm,b- Because this will be used (see end of Section 4.6) 
more for aesthetic than essential reasons, we take the opportunity to omit 
the proof. 

Lemma 15. cm+i,b <cm,b +{M + 1)G'{c*,B), where G' {c* , B) depends 
only on c* and B. 

3.3. The upper hound. Proposition 17 below gives the upper bound in 
Theorem 2. If the flows in Proposition 5 were exactly standardized global 
flows, then Proposition 17 would follow immediately. Instead, we will use 
the following easy smoothing lemma, together with the patching lemmas of 
Section 3.2. 

Lemma 16. Let take values in the N x N torus. Let Ul be uniform 
on the Lx L square centered at the origin (L odd), independent of Z^. Sup- 
pose N~^Zj\f converges in distribution to the uniform probability distribution 
on [0,1)^. Then, provided Lj\f /N ^ sufficiently slowly, 

min P{Zn + Ul^ = v)/N-^ ^ 1. 

Proposition 17. There exist standardized global flows f(-|c) on the N x 
N torus such that maxe/(e|c) < B and such that, for each fixed M, 

limsupAf~^EcostAr(f(-|c),c) < M~'^cm b- 
N 
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Proof. Fix M\N and an environment c on the N x N torus. We con- 
struct a flow in 4 steps, where only step 2 depends on c. Consider the 
transportation measure pn,m from Proposition 5. The restriction of its en- 
trance marginal to the boundary Bouj^^) of an extended M x M square S in 
the natural partition is (up to translation) a measure 9 on Boum, which (by 
Properties 3) does not depend on S and is the same for the exit marginal. 

Step 1. Each vertex v of the N x N torus is in some M x M square S 
in the natural partition. Construct a path- flow (with some flow- volume fi) 
from each f to a random position (bo, say) on the boundary Bou(5') with 
distribution proportional to 6, with total volume A'^"^ starting from each v. 

Step 2. Use the path-flow pn,m{'\c) in Proposition 5 [with flow-volume 
f2 = f2('|c)] to send flow from a typical position 6o ^ Bou7v,a/ to another 
position bi € Bouat^m- 

Step 3. Reverse step 1, sending flow from bi to a uniform position {v' , 
say) inside an adjacent M x M square. Write fs for the flow-volume. 

Step 4. Send flow from v' to a uniform random position in the L^v x L^v 
square centered at v' . Write £4 for the flow- volume. 

Let f = fi -|- f2 + fa + f4 be the flow- volume and H the transportation mea- 
sure for the concatenated flow. Combining Proposition 5(iv) and Lemma 16, 
provided Ljv/A^ — > slowly, we have 

(26) min E{yo, yi) = N^^wn, 

yo,yi&T^ 

where ■wat— >lasA^— >oo. In step 1, in each square 5 we are in the setting of 
Lemma 8, seeking a path- flow whose transportation measure is of the form 
p X p, where p has total mass N~^M'^ . By Lemma 8, we can choose such a 
flow to satisfy 

max/i(e) < 2maxp(6) < 2iV~U'/^ 

e b 

The same bound holds for £3 in step 3. In step 4, using Lemma 7 and scaling, 
we can take 

max /4(e) <iV.^ = —. 

So we have 

(ff^^ff^^ff 4M2-K1/4)L^ 
max(/i(e) + /3(e) + /4(e)) < — , 

while, by Proposition 5, 

max /2(e) < B. 

e 

We want to use Lemma 13. Fix 6 < B — j. Recall (26): by deleting flow, we 
may assume 



FLOW THROUGH DISORDERED LATTICE 25 

So w^^f is a standardized global flow. Apply Lemma 13 to w'^^f = w]^^i2 + 
t/;^^(fi + fs + U) with i?o = w]^B. For sufficiently large N , we have 

<d; d< -; b<B 

and so Lemma 13 gives a standardized global flow f (-Ic) such that maxe /(e|c) < B 
and 

-2 / , ^ N 85N^B^w7r^ ^ 

COSt(Ar)(f(-|c),c) < COSt(jv)(f2(-|c),c) + ^^^^ _ . 

So using Proposition 5(ii), 

N ]Ecost(^)(f(.|c),c)<^ + ^^— 

Since wn ^ ^ and (5 > is arbitrary, we see 

limsupiV-2EcostAr(f(-|c),c) < M-'^cm b 

N 

as N runs through multiples of M. Using Lemma 14 to interpolate other 
values of iV, we have established Proposition 17. □ 

4. The lower bound. In outline, we will use a concentration inequality 
(28) and a convexity property (Lemma 18) to derive a rather technical lower 
bound (Lemma 19) on the cost of any flow across a M x M square in terms of 
a smoothed transportation measure across a slightly larger square. We then 
consider properties of the induced flow across subsquares arising from a 
standardized global flow (Section 4.4). These ingredients are then combined 
to prove the lower bound of Proposition 22. 

4.1. A concentration inequality. Return to the setting (Section 2.4) of 
the extended M x M square. Consider environments c satisfying 

(27) < c(e) < c* Ve. 

Given a transportation measure Q on Bouj\/ x Bouj\/, and an environment 
c, we defined 

costM,B{Q, c) = inf < costM(flo(^), c) : tra(/i) = Q, maxflo(^)(e) < B 

By compactness, the inf is attained. Now consider the effect of changing c 
to c by changing the cost-factor on only a single edge cq. If ^ attains the 
inf above, then 

costAf (flo(/_f), c) < costM(fio(^),c) + c* B^ 
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and it follows that 



COStA/,B(<5,c) - COStM,BiQ,c) \ < C* B"^ . 



So in our model (3) of i.i.d. random environment we can apply the most basic 
form of the method of hounded differences (Azuma-Hoeffding inequality: see, 
e.g., [18], Section 1.3 or [5], Section 7.2) to conclude that, for A > 0, 



Note that \£m \ = 2M^ [recall convention below definition (10)]. Note also 
that we only need a one-sided bound. 

Discussion. The expectation in (28) is order M^, and the inequality says 
that fluctuations are only of order Af . If we use the same approach to the 
kind of "maximal flow subject to i.i.d. edge-capacities" problem mentioned 
in the Introduction, then we would study a r.v. measuring maximal flow 
across a M x M square. As with (28), this r.v. would have fluctuations of 
order M, but now the expectation is also of order M, so we do not have an 
easy concentration result. 

Inequality (28) is also the reason we impose the bound B on edge-flows. 
This bound causes complications in Section 4.3, where to "smooth" trans- 
portation measures on the M x M square we are forced (to preserve the 
bound B) to extend to a larger square. 

4.2. A convexity lemma. For fixed c the function (5) f costA/(f, c), is 
convex, so we have the following: 

Lemma 18. For fixed c, the map Q — > costM,B(Q, c) is convex, as a 
function of feasible transportation measures Q on Boujv/ x Boujv/- 

4.3. Smoothing the transportation measure. Inequality (28) refers to a 
fixed transportation measure Q. We need to use it when Q depends on 
c. We do this by defining a certain finite set [Ss{M + 2K, K, B) below] 
of transportation measures and then comparing a general Q to a nearby 
element of that set. 

Fix K\M. Consider the extended square [0,M — 1]^ centered inside the 
larger extended square [—K,M + K— 1]^. Let /i be a path-flow across [0,M — 
1]^ which is feasible in the sense 



so that Q := tra(^) is a measure on Boum x Boua/ . We will first use Lemma 10 
to construct an extension to a path-flow fiext on [—K,M + K — 1]'^ whose 
transportation measure Qcxt := ti^a(/icxt) has a smoothness property. 



(28) P(costM,B(Q,c) < -A + EcostM,s(Q,c 




maxflo(/i)(e) < B 
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Fig. 3. Boundary points of an extended M x M square within matching boundary points 
of an extended (M + 2K) x (M + 2K) square. Here M = 1Q,K = 2. A typical K x K 
square is shown at the top. The lines on the right side indicate part of the partition of 
those boundary vertices into blocks of size K. 

Remark. As indicated in Section 2.5, it would be more natural and 
elegant to do the construction of this section by rerouting flow within the 
M X M square. But we are unable to do so. The use of the larger square 
[— iC, M + ET — 1]^ is purely a technical device, allowing us to create flows with 
a smooth transportation measure on the larger square without rerouting flow 
in the inner square. 

The construction is illustrated in Figure 3. Partition the 4M boundary 
points into {AM)/K adjacent blocks of size writing block(6) for an index 
of the block containing h. Consider a K x K square R which extends from 
a block in the boundary of the M x M square to the parallel block in the 
boundary of the (M + 2K) x (M + 2K) square (see Figure 3). Write fimner 
and i?outcr for the inner and outer boundaries of R. Write 9i and 02 for the 
restrictions of Qcnt and of Qcxi to i?inner- By feasibility, ^i(6) +02ib) < B for 
b € -Rinner- Lemma 10 shows that we can create a feasible flow through R 
whose transportation measure is 62 x 02 + 61 x 9i, where 0i is the constant 
measure on -Router whose total mass equals ^i(i?mner)- Doing this for each 
square R constructs a feasible flow /Xsmooth across [-K,M + K -1]'^. Because 
each square R has 2K'^ edges, and there are AM/K squares, and each edge 
cost is at most c*B^, we have 

COStAf+2i^(flo(^smooth), C) < COSt A/ (flo(;U) , c) + 8c*B'^MK. 

We emphasize this holds for any path-flow /z and any environment c satisfy- 
ing (27); so the flow may depend on the environment. Note that Qsmooth '■= 
tra(^smooth) depends only on Q :=tra(/z). Recalling the definition (9) of 
costM,B(QiC), we can minimize over path-flows with prescribed transporta- 
tion measures to get 

(29) COStA.f+2if,B(Qsmooth,c) < COStAf,B(Q, c) +8c*B'^MK. 
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Recall block(6) denotes the block containing b. Define the set S{M + 2K, K, B) 
to consist of transportation measures Q on Bouj\/+2/< x Boujv/+2_ft' with the 
properties: 

(i) Q{b,b') depends only on block(6) and block(6'). 

(ii) Q is a feasible transportation measure for some path-flow fx across 
[-K,M + K -1]^ with maXeflo(/i)(e) < B. 

(iii) Q{b,b') = if either b or b' is not one of the 4M boundary points 
which are parallel to the 4M boundary points of the M x M square. 

By construction, Qsmooth := tra(//smooth) is an element of S{M+2K, K, B). 
We want to relate drift (Qsmooth) to di\ii{Q). Write R for a typical K x K 
square, Bou|/j for its boundary and Q^ji for the transportation measure 
across R, in the construction above. Then 

E (b' -b)QsmootUb,b')= E {b' -b)Q{b,b') 

b,b' GBovlm+2K bjfe'GBoujv/ 

+ E E {b'-b)Q\nib,b'). 

R 6,fe'eBou|fl 

For fixed R, writing and for the total Qcxi-measure and the Qent- 
measure of the inner boundary of R (the co-boundary with the M x M 
square), we have (using the product form of (5|_r) 

E ib'-b)Q^R{b,b') 

6,6'GBou|fl 

< \q^-q^\K+{q^ + q^){K/2) 
<lK{q^ + q^). 

Summing over the squares i?, the size of the second term on the right in 
(30) is at most Miq, where q:= Q(Bouj\,/ x Boua/). So 

||(M + 2K)2drift(Qsmooth) -M2drift(Q)||i <3Kg. 
Noting that || drift(Q)||i < |f and that 1 - ,^M+2Ky^ - '^K/M, we find 

II dntt(Q_„) - dntt(Q)||, < + ^^jj^ < 

Note that, by feasibility, 

q < 2MB. 

Now let (5 > and define SsiM + 2K,K,B) to be the subset of S{M + 
2K, K, B) consisting of those Q which satisfy also 

(iv) Q{b, b') is a multiple oi 5/K. 
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We first show that the cardinahty of Ss{M + 2K, K, B) is bounded by a 
number depending only on M/K and B /5: 

(31) \Ss{M + 2K,K,B)\<l{M/K,B/6):=[l + j\ 

To see this, remove constraint (ii) from the definition of S{M + 2K,K,B). 
Then a Q satisfying (i), (iii), (iv) is determined by an array (m(/3i,/32)) of 
nonnegative integers, where (3i and (32 run through the AM/K blocks and 

Q(6i,62) = m(/3i,/32)W he Pi. 

The latter implies that the volume of flow from bi to /J2 is m{f3\,f32)6, and 
because of the bound B on edge-flows, we must have m(/3i,/32) < B/5. This 
establishes (31). 

Define a "trimmed" transportation measure Qtrim < Qsmooth by 
Qtrim(^5&') is the largest multiple of 5 smaller than (3smooth(&i ^')- 

So Qtrim is an element of Ss{M + 2K, K, B). 
Now consider a path-flow v with 

tra(z^) = Qtrim! maxflo(z^)(e) < B. 

e 

Note that Q := Qsmooth — Qtrim has Q(&i,&2) < ^/K and so its marginals 
satisfy max(Qcnt(^)> '3cxi(^)) < AM5/K. Using Lemma 9, we can construct 
a path-flow v' such that 

tra(z^') = Q; maxflo(z/')(e) < 16M(5/ii: 

and so 

tra(z/ + I/') = Qsmooth; flo(z^ + J^')(e) < flo(i^)(e) + IQM5/K < B' , 

where we set B' := B -\- 16M6/K. Because increasing the flow- volume across 
an edge by r] (subject to the bound B') can increase the cost associated with 
that edge by at most 2c*B'r], it follows that 

costM+2ft'(flo(z^ + u'),c) < costA./+2ft'(flo(i^),c) + 2B' c* {IQM5 / K)\£m+2k\- 
Minimizing over path-flows with prescribed transportation measures, we see 

COSiM+2K,B' (Qsmooth, c) 

(32) 

< cost A/+2i^,B (Qtrim, c) + UB'c*6M{M + 2K f/K. 
We now bring into play the probability measure on environments. Write 
(33) mB((5) :=EcostM+2i^,B(Q,c). 
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Consider, for each c, some feasible path-flow /i(-|c) across the extended 
M X M square (this notation emphasizes we allow dependence on c); write 
Qsmooth('|c) and (5trim('|c) for its smoothed and trimmed transportation 
measures. From the concentration inequality (28) applied to M + 2K, and 
the finiteness of Ss{M + 2K, K, B), 

P(cOStM+2i^,B(Qtrim(-|c),c) < - A + ((5trim(- |c))) 

<WA-,B/i)exp(-^5^^^^^^^)^ 

This remains true if we replace the first occurrence of Qtrim by the larger 
Qsmooth- Combining with (29) and rescaling A, 

F(costM,B(Q(-|c), c) < -(8 + \)c*B^MK + mB{Qtrim{-\c))) 

<f(M/A-,BM)exp(-^^^j^). 

We tidy slightly by setting K = M/L for some L > 4 and setting A = 4 to 
get 

P(costM,B(Q(-|c),c) < -120*3"^ M^/L + mB{QtrMc))) 
<£{L,B/6)exp{-M^/L'^). 

Note that for any feasible Q there is an upper bound on niBiQ) implied by 
feasibility, which works out as mB{Qtnm{-\c)) < B^c*2{M + 2K)'^ . Applying 
the elementary fact, for nonnegative r.v.'s X,Y, 

if P{X < y - a) < 6 and y < 2/0, then EX>EY -a- byo 

shows 

EcostA./,ij(Q(-|c),c) >EmB(Qtrim(-|c)) - 12c*B^M^/L 

- B'^c*2{M + 2Kf£{L,B/d) exp{-M'^ / L'^). 

Now (32) says 

mB{Qtnm{-\c),c)>mB'{Qsmooth{-\c),c)-64B'c*6M{M + 2Kf/K. 
Convexity (Lemma 18) and Jensen's inequality imply 

EmB'(<5smooth(-|c)) > "lB'(E(5smooth(-|c)). 

Combining the last three inequalities and using the fact M + 2K < 
(because L > 4), 

EcostM,B(Q(-|c),c) > mB,(IEQsmooth(-|c)) - I2c* B^ — 

- lUB'c*SM^L - bcB^M^£{L,B/S)exp{-M^/L'^). 
We can summarize the results of this section as the next lemma: 
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Lemma 19. Let M = KL, where L > 4. For each environment c with 
maxec(e) < c* , let At(-|c) he a path-flow across the extended M x M square 
with maXeflo(/i(-|c))(e) < B. Let Q{-\c) be its transportation measure and 
Qsmooth('|c) the smoothed transportation measure across the extended (M + 
2K) X (M + 2L<:) square. Then 

(34) ||drift(Q,„ooth(-|c)) -drift(Q(-|c))||i < 
and for any 5 > 0, 

(35) > M-2mB'(EQsmooth(-|c)) - 12c* B^L 

- U4B'c*6L - 5c*B^£{L,B/6) exp{-M^ /L'^), 
where B' = B + 166L. 

4.4. Local flows induced by a global flow. This section derives identities 
(Lemma 20) which restate certain quantities for flows on the N x N torus in 
terms of induced empirical distributions across M x M squares. Of course, 
this is "just bookkeeping" ; the point is to obtain quantities defined on spaces 
depending on M not on N, so that it makes sense later to let N oo. Until 
further notice, fix 2 < M < A^, a global environment c on the N x N torus 
and a standardized global flow fj,. We start with some observations. 

Observation 1. A path vr in the N x N torus from to induces 
(by using the same increments) a path in the lattice from to a point 
z* of the form z* = z'^ -\- {iN,jN), where the integer pair {i,j) indicates net 
winding around the torus. Define 

(36) iON{Tr) = {z* - z^)/N eR^ 

so that u;Ar(7r) indicates both the normalized relative positions of origin and 
destination of the path within the torus, and also the number of windings. 

Observation 2. Write Ai{S) for the space of measures on a space S. 
If -0 is a probability distribution on A4{S), then (under the natural integra- 
bility condition) it has a mean measure ■0 G A4{S): 



Jm{s) 

Observation 3. Given a measure v ^ M.{S x S') with first marginal 
vi G Al(5), there is for s G 5 a conditional (probability) distribution i'{-\s) G 
M{S') satisfying 



u{A xB)= I v{B\s)ui[ds). 
J A 
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Recall that Sm and Cm denote the sets of paths a and of environ- 
ments c on the extended M x M square. Given (x, y) € [0, — 1]^, consider 
[x,x + M — 1] X [y,y + M — 1] as a subsquare of the N x N torus. The 
global environment c induces a local environment on the subsquare which, 
by translation, becomes an element c^^y = {cx,y{e)) of Cm- 

(37) Cx,yie) =c{e + {x,y)), cGSm, 

where e + {x,y) denotes the edge e translated by {x,y). 

Consider a path vr in the N x N torus. Then vr intersects the subsquare 
[x,x + M — 1] X [y,y + M — 1] via some number (maybe zero) of segments 
a, which (by translation) can be regarded as elements of T,m- Call these 
segments o"i(7r; x, y), say. So the standardized global flow n induces a measure 
fj,x,y on Sa/ via 

^^x,y{■) =^^^{'^■cri{^T;x,y) G •}. 

i 

We need to record also the value of u!i\[{tt) from (36), which leads us to define 
the measure fi^y on Sm x M^: 

fJ'tyi- X ■)=Y^fJ'{TT:ai{TT;x,y) G •,u;7v(7r) G •}. 

i 

Note 

{cx,y, fJ-t,y) is an element of Cj\/ x A^(Sa/ x M^). 

Now take (x, y) as a uniform random vertex of the N x N torus, so that the 
random element {cx,y, IJ^t y) has some probability distribution, which we call 
r. So 

(38) r is a probability measure on Cm x J^i^M x I^^)- 

In words, T describes the empirical distribution of local (i.e., across M x M 
squares) flows, jointly with the local environment and the relative source- 
destination vector of flow-paths. As indicated initially, F contains quite a lot 
of information about the global flow and environment, /i and c, as we next 
describe. 

The identity for costs. Write i^^ ^ v for the "take first marginal" map 
7W(Sa/ X ]R2) ^ M{^m)- Then 

(C, U^) COStA/(flo(z^),c) 

is a functional defined on Cm x A4(T,m x M^). We assert 

Af~^ / costA/(flo(z^),c)r(dc,(ii/+) =iV~^cost(flo(^),c). 

JCmxM{SmxR'^) 

(39) 
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This holds because the integral equals 

N^'^ ^ COStM(flo(/ia:,y),Ca:,y) 

and the sum here equals cost(flo(^), c) because for each edge e G £[n)i 
the cost c(e)/^(e) is counted in exactly subsquares. 

The identity for drift. As in Observation 1, consider a path vr in the 
discrete torus [0, A'^ — 1]^ from to and the induced path in the lattice 

from to a point z* . Consider the natural partition IIo of [0, — 1]^ 
into iV^/M^ subsquares S^^y = [x, x + M - 1] x [y, + M - 1] . Then 

where the second sum is over the path fragments a = ai{7r;x,y), where vr 
intersects Sx^y Averaging over shifts of the partition gives 

(40) z*-z^=M-^ J2 EM^)-ent(cj)). 

{x,y)e[0,N~l]'^ " 

Now let r2 be the second marginal of T. So r2 is a probability distribution 
on A^(Sa/ X M^). Let r2 be its mean measure (Observation 2). Concretely, 

^2 = N~^Y.x,y^4,y For ZGZ2, 



Using linearity of the map u — > drift(tra(z^)), 

z ' 



drift! tral Tal • x ^ 



(41) 



M-2Ar-2^ /'^(exi(a,(^;x,y)) 

x,y i 

-ent(cri(7r;x,y)))l(^^(^)=2/jv)/i((i7r) 
A^"^z//|7r:u;Ar(7r) = 



the last equality by (40). Write 

T,v(-)=^^~V{vr:cuyv(vr)G-}. 
From the definition of standardized global flow, 

the push- forward of T^v under the map u^u mod (1,1) 

(42) 

is the uniform probability distribution on{A^ "'^z:zG[0,A^ — 1]^}. 
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Note that (41) can be rewritten as 

(43) drift(tra(f 2(- x {u}))) = mT7v(u), u G R^. 

Identity for marginals of transportation measure. When a segment (Tj(7r; x, y) 
exits the square [x, x + M — 1] x [y,y + M — 1] via some boundary vertex, it 
must enter an adjacent square. This fact easily imphes the following identity. 
Given u = z/N G R^^ ^j-j^-g q ^ tra(r2(- x {u})). Then 

the push-forward of Qcxi (restricted to Boum) 

(44) under the map h ^reflect 
equals Qent (restricted to Boum)- 

Identity for terminal vertices. For a path a € Sm, let inteM('7) G {0, 1, 2} 
be the number of endpoints [ent(c7) and exi((T)] which are in the interior 
B'j^ rather than the boundary Boua/ . Here inte is a mnemonic for interior. 
Because each path in S^jy) has exactly two endpoints, and an endpoint 
appears in exactly subsquares, 

^ I inteA/(cT)/x^,j,(da) = 2MV(S(^)) = 2M'^N. 

x,y -^^jv/ 

Writing f2i(-) = N~'^J2x,y t^x,y{-), so that is the marginal of r2 above, 
we find 

(45) / mteM{(T)r2i{da)=2M^/N. 

Bound for windings. For a path vr in the N x N torus, 

n(7r,e)>iV||a;;v(vr)||i 

and so for a standardized global flow /x with flow-volume f , 
E f{e)>N j ||w7v(vr)||iMd^). 

If the capacity constraint maxe/(e) <B is satisfied, then the left-hand side 
is at most 2N^B and so 

(46) J ||u;Ar(7r)||i fj,{dir) < 2NB. 

Identity for total mass. From the definitions we can write the total mass 
of r2 as 

f2($]M xM2)=Ar-2^ f H.,,y{7T)^x{d7T), 
x,y 
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where Hx^y{ir) is the number of segments of vr which intersect the square 
[x,x -\- M — 1] X [y,y + M — 1]. Because each step enters exactly M such 
squares, and ent(7r) is in exactly such squares, 

(47) f 2(Sa/ X M^) = J^M^iV + M J len(7r) fiidn)^ , 

where len(7r) is the length of the path vr and where the factor N is the total 
mass of jj,. 

So far we have worked with a fixed standardized global flow fi and a 
fixed global environment c, so rewrite the T at (38) as T^'^ to indicate 
this explicitly. Now suppose the environment c is random according to our 
probability model (3), let /^(-Ic) be a standardized global flow depending on 
the realization of c, and write 

(48) r(-)=Er^(-i^)'^(-) 

for the mixed probability distribution on Cm x A1(Sj\/ x R^). Note that the 
marginal (Fi, say) on Cm is clearly just the i.i.d. distribution (3). Note also 
that (39) implies 

M"^ [ costA/(flo(i/),c)r(dc,dz^+) 

JCm X (Em xIR2) 

(49) 

= iV"2Ecost(flo(/i(-|c)),c). 

We now want to replace the probability measure F by a (nonprobability) 
measure f on Cm ^ ^ defined (loosely speaking) by replacing the 
conditional distribution (given environment) on path-flows by its conditional 
mean. Precisely, for c G Cm, the probability distribution T has an associated 
conditional distribution Tc on M(T,m x M^). Define fc G MiT,^ x M^) as 
the mean measure (Observation 2) of Fc, and then define 

(50) f(dc,-,-)=ri(dc)fc(-,-)- 

Lemma 20 gives identities and inequalities relating f to the underlying global 
flows in the random environment. Identities (43)-(45) and bound (46), which 
do not involve interaction between flows and environment, extend to the 
random environment setting in obvious ways, recorded as (iii)-(vi) below. 
We handle (49) as follows. The map 

(51) ^ 1/ ^ flo(z^) from M(J:m x^'^)^M(J:m)^:Fm 

takes fc to some element of Tm, say, fc. Applying Jensen's inequality to the 
convex function — > costA/(flo(z^), c) and the probability measure Fc gives 

costjv/(fc,c) < / costM(flo(z^),c) Fc(dz^'^). 

JA4{SmxR2) 

Integrating both sides against Ti{dc) and using (49) gives the inequality 
stated in (ii) below. 
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Lemma 20. For arbitrary standardized global flows /i(-|c) on the N x N 
torus, define the measure f on Cm x x by (48), (50); 

(i) T{dc,-,-) =Hi((ic) rc(-,-), where Hi is the i.i.d. probability distri- 
bution (3) on Cm cLn-d where G Ai{T,M x 1^^) for each c € Cm- 

(ii) For fc defined below (51), 

M-2EcostA/(fc,c) < A^"^Ecost(flo(^(-|c)),c). 

(iii) Write f 2(-, •) = f (Cjv/ x • x •) for the marginal measure on Sjv/ x M^. 
For I? C wnie Qt-^l = tra(f 2(- x D)). Then 

drift(Q[^])= / uTN{du), 

JD 

where 

Tjv(-) =^^~^E^{v^:a;^(^) S -Ic}. 

(iv) ForDcM.'^, 

the push-forward of Q^^} (restricted to Boum) 
under the map b —>■ b^^^'^'^^ 

equals (restricted to Boujv/J. 

(v) Forr2i(-) = ri(- x]R2), 

f inteM(fT) f 2i((i(T) = 2M^/N. 
For the remaining parts assume maxegf^j^^ flo(/u(-|c))(e) < B\/c. 

(vi) /iR2 ||n||iT7v(dn) < 2B. 

(vii) maxeef^, /c(e) < BVc. 

(viii) r(CA/ X Em x M^) < Af2^-i + 2MB. 

Proof. Part (i) holds by definition and parts (ii)-(vi) were discussed 
above. Part (vii) is straightforward. For part (viii), by (47), 

f (Cm X Sa/ X M^) = M^AT-i + N'''^ME J len(7r) /x(d7r|c). 

Now the integral equals 

J2 flo(M-|£))(e)<S|f(7V)l, 

ee£'(]v) 

giving (viii). □ 
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4.5. Properties of weak limits. Lemma 20 dealt with arbitrary standard- 
ized global flows /u(-|c). So we can apply it to flows attaining the inf in the 
deflnition 

cost(^-)(c, := inf{cost(f , c) :f a standardized global flow satisfying (4)} 

and then see what happens in the N ^ oo limit. 
Define 

(52) 7*iK,B) :=linynfiV-2Ecost(^)(c,B). 

Let C be the subset of paths a which go across the M x M square 
in the sense that both ent{a) and exi(o") are in the boundary Boum- 

Lemma 21. Fix M >2. There exist measures Hc(-, •), c G Cm on x 

M? with the following properties. Let fc be the flow-volume for the path-flow 
Hc(- xM2).. 

(53) max/c(e)<S Vc, 

(54) A^r'^E cost Af (fc , c) < 7* (k, 5) . 

For bounded D C M^ write Q^^^ = tra(EHc(- x D)). Then 

(55) e Qm, 

(56) drift(Q[^l)= / uT{du), 

Jd 

where T is a probability measure on M? whose push-forward under the map 
u mod (1, 1) is the uniform probability distribution on T^. 

Note that the final sentence implies (via a Radon-Nikodym argument) 
that Q :=tra(EHc(- x M^)) has a "density representation" 

Q = I Q["lT((it/), 

where 

QMsQm; drift(QM) = n, 

which implies Q € Qjv/ • This result is the fundamental motivation for con- 
sidering Qm ■ 

Proof of Lemma 21. RecaU the vague topology on M{Cm x x 
M^), that is, convergence of integrals of continuous functions with compact 
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support. Recall that in this topology the condition for relative compactness 
of (Tn) is 

sup Fat (J^) < oo yK compact. 

N 

Recall also the Fatou-like result: 

if r^r — > H vaguely, then J hd'B< lirninf J hdTjy 

for continuous h>0. 

For each A^, we consider the standardized global flows attaining the inf in the 
definition of cost(7v) (c, B) and then define r^v to be the measure T associated 
with these flows as in Lemma 20. The relative compactness condition holds 
by (viii) of Lemma 20, so we can flrst take a subsequence of N where the 
lim inf in (52) obtains, and then a subsequence such that Ftv ^ S weakly, 
for some measure H on 

Cm X Sm X ]R2. By (v) of Lemma 20, 

J inteA/((T) H(dc, da, du) = 0, 

which says that H is supported on Cm x x M^, that is, we may replace 
T,M by S^. For each N, the marginal of F^v on Cm is the i.i.d. law Hi, so 
the marginal of H is the same law Hi, and so there is a representation 

E{dc,-,-) =Ei{dc)Ec{-,-) 

which serves to define Hc(-,-). Properties (53)-(56) now follow as ^ oo 
limits of the corresponding properties (ii), (iii), (iv), (vii) of Lemma 20, once 
we have shown that the probability measure Tjv on converges weakly to 
some limit T. But (vi) of Lemma 20 implies tightness, so by passing to a 
further subsequence in A'", we may assume Tat — > some T. □ 

4.6. Completing the lower bound. Now we will combine previous ingre- 
dients to show the following: 

Proposition 22. For any B" > B, 

lim sup M~^cm _B" <7*(k,-B) := liminf A^^^Ecost(7v)(c, -B). 

M ' ^ 

Lemma 11 shows 

B "f*{K, B) is continuous nonincreasing on ^<B < oo. 

Now Propositions 17 and 22 combine to prove Theorem 2 and, hence, The- 
orem 1. 
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Proof of Proposition 22. The proof mostly consists of combining 
Lemmas 19 and 21. Fix 

M = KL; 5>0; B>1 

and define 

e = 22B/L; B' = B + 16SL 

and let B" >B' + e. First fix D C and consider 

Q[^1(-|c) :=tra(He(- xI?)), 

where Hc(-, •) is given in Lemma 21. Applying the linear map Q — > Qsmooth 

to Q^^^{-\c) gives transportation measures Qlmooth('l'^) across the extended 
(M + 2K) X (M + 2K) square such that 

d"ft(Ql2ooth(-|c)) - drift(Q[^l(.|c)) = {x,,D,yc,D), 

where, by (34), |xc,d| + I^c.dI < £• If (say) Xc^d > and yc,D < 0, then in the 
notation of Lemma 12 (applied to M + 2K) the "adjusted" transportation 
measure 

«3Ldj(-|c) :=doth(-|c)+^c,DQ^ + yc,DQ^ 

has 

(57) drift(Q2j(.|c)) = drift(Q[^](.|c)). 

The effect of the general Q Qsmooth map on marginals is simply to average 
over blocks, and so the fact (55) 

q[^1:=EQ[^1(.|c) gQm 

implies 

Qlmloth -JEQSooth(-lc) e Qm+2K- 

So now 

<?!sj-E<3a(-ic) 

is of the form 

where J2i£i,D < £■ So clearly Qj^^j G Qm+2K- And 
drift(Q2j)=Edrift(Q2j(-|c)) 



Edrift(Q^^H-lc)) by (57) 
drift(Q[^]) 

/ uT{du) by (56). 

JD 
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Writing Qadj for the D = M.'^ case of Q^J ) this density representation (cf. note 
below Lemma 21) shows Qadj £ Qm+2K- ^PPly Lemma 12 with (Qadj i Qsmooth := 
Qlmooth, B', M + 2K) in place of (Q', Q, B\ B, M): 
(58) EcOStAf+2i^,B'(Qsmooth, c) > Cm+2K,B" " Ac* B" e{M + 2Kf. 
Recalling the definition (33) of mB{Q), 

TUB' (EQsmooth ( ■ |c) ) = IE COStM+2K,B' (Qsmooth (" I c) , c) 
= E COStM+2i^,_B' (Qsmooth , c) . 

Now we apply Lemma 19 to the measures Hc(- x •) given by Lemma 21, with 
fc the associated flow-volume; conclusion (35) becomes 

M-2EcOStM(fc,c) > M-2EcOStM+2J^,B'(Qsmooth,c) - 12c* B^ / L 



The left-hand side is upper bounded by j^{k,B) by (54). So combining this 
inequality with (58) gives 



the conclusion of Proposition 22, except with the limsupj;// taken through 
a subsequence Mj with Mj^i/Mj — > 1. But then Lemma 15 identifies the 
limsup through that subsequence with the actual limsup. □ 

Remark. Without using Lemma 15, the argument gives the analog of 
Proposition 22 with liminfA// in place of limsup^, which in turn gives the 
(less aesthetic) analog of Theorem 2 with liminfM in place of lim^- 

5. Discussion. We start with comments tied closely to the statement and 
proof of Theorem 1 and then venture further afield. 

1. We can apply the "method of bounded differences" argument (28) 
directly to the random cost cost(jv)(c, S) in Theorem 1, and conclude 



This of course implies the SLLN corresponding to Theorem 1. Various other 
"optimization over random data" problems have this character — that it is 
easy to show concentration about the mean, but not easy to show existence 
of a limit constant for the mean. Perhaps best known is random 3-SAT [8]. 



U4B'c*6L - 5c*B^e{L,B/6) exp{-M^/L^). 




P(|cost(Af)(c,S) -Ecost(Ar)(c,5)| > A) < 2 exp 
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2. Our setup for Theorem 1 did not exclude the possibihty P(c(e) = 
0) > 0. Recall [6, 12] that the critical value for bond percolation in 1? equals 
1/2. It seems intuitively clear that 

if P(c(e) =0) > i, then 7(k,S) = 0; 
if P(c(e) =0) < i, then 7(k,S) >0, 

but we have not checked carefully. 

3. The limit 7(k, -B) is a priori nonincr easing in B. By considering ran- 
dom minimum-length paths (Lemma 6) and the bound c* on the support of 
K, we find an upper bound 

7(k,5) <c78. 

It is not hard to formalize the intuitive idea that as | 1/4 the only feasible 
flows have /(e) ~ 1/4 for almost all e, and so 

lim 7(k,5) = iEc(e). 

4. Relaxing the requirement of bounded support, for what k do we still 
expect a finite limit ^{k,B)? This question seems rather subtle. It is not 
obvious that the condition Ec(e) < oo is necessary, because we might be 
able to route flow to avoid edges e with large c(e). On the other hand, the 
condition Ec^/^(e) < cxo is necessary, otherwise [by considering the minimum 
of c(e) over the four edges at a vertex] we have Ecost(7v)(c, B) = oo for finite 
N. 

5. Return to the case where k has bounded support. Define cost^jy) (c, oo) 
as the "-B = oo" case where there is no bound on edge capacity. By mono- 
tonicity (of limit constants in B), the limit constant 

7(k,oo):= lim j{k,B) 

B—iOO 

exists. From Theorem 1 we get a bound 

limsupA^~^Ecost(^-)(c, oo) < 7(^,00). 

N 

Conjecture 23. limAriV~^Ecost(Ar)(c, 00) =7(^,00) = limM->oo x 

CM, 00 ■ 

It seems plausible that one could just modify the arguments in this paper 
to directly study the B = 00 case. Loosely speaking, in the B = 00 case we 
expect that the optimal flow /(e) across a typical edge e has exponential 
tail 

(59) P(/(e) >x) <^iexp(-2;/^2), < a; < 00, 

so that imposing a large bound B has little effect. 
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6. It is tempting to speculate that the Q € Qj^f which attains the inf 
in the definition (10) should be some "geometrically natural" measure not 
depending on k or S. One possibility: take the intersection of Euclidean- 
isotropic lines with the continuous square [0, M]^ to get a transportation 
measure on the boundary of this square; then discretize to get some isotropic 
transportation measure Qm on the boundary of the discrete M x M square. 

7. The construction in Section 3 gives flows of asymptotically optimal 
cost. From the use of the law of large numbers in the proof of Proposition 4, 
one can see that the routes of this flow (after scaling the N x N torus to the 
continuous torus [0, 1)^) converge to straight line segments in the continuous 
torus. But we have not proved these are the minimum-length straight line 
segments; we have not excluded the possibility that some non- vanishing 
proportion of flow takes a long route around the torus instead of taking the 
shortest route. This possibility in taken into account in Lemma 21 (and the 
definition of Qf/), where we use a general T rather than just the uniform 
distribution on [-1/2, 1/2]^. 

8. Given an environment c on the N x N torus, there is a calculus- 
type condition for a standardized global flow f to be a local minimum of 
the function f ^ cost(jv)(f,c). This condition (Waldrop equilibrium [11]) is: 
for each source-destination pair (x,y), the flow from x to y can only use 
minimum-weight paths, where here the weight of an edge e is c(e)/(e). But 
we do not know how to exploit this condition in studying the optimal flow. 

9. If instead of the quadratic costs (1) we used linear costs Z)eC(e)/(e), 
then the optimal flow simply chooses the minimum-cost path [where cost of 
a path equals sum of edge-costs c(e)] between each source and each destina- 
tion; taking B = oo, there is no interaction between flows. This case is just 
first passage percolation, and mean costs in our "flow" setting are easily re- 
lated to the time constant in first passage percolation. Questions concerning 
flow volume across edges have apparently not been studied rigorously, but 
one expects power-law tail behavior instead of (59). 

10. There is a huge literature on algorithms for flows in networks, illus- 
trated by the monograph [1]. But the kind of multicommodity flow prob- 
lems that we study are typically NP-hard as algorithmic problems. There is 
a large body of theoretical work going back to Leighton-Rao [15] relating 
multicommodity flow to other network problems, but this typically gives up- 
per and lower bounds differing by log(number of vertices) factors. We should 
emphasize that the basic max-flow min-cut theorem does not apply in our 
multicommodity flow setting. 

11. The kind of feasible flow problems mentioned in the introduction 
have been studied in the unidirectional case where one seeks to maximize 
volume of flow from somewhere on the bottom side of a square to somewhere 
on the top side. See [7, 9, 13, 21] for results and connections with percolation 
and first-passage percolation. In this setting one can apply the max-fiow 
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mill-cut theorem. Otherwise there seems no Hterature very closely related 
to Theorem 1, though of course the general idea "consider large blocks" is 
pervasive throughout the study of spatial stochastic processes. 

12. A more abstract view of our argument goes as follows. Subsequential 
weak limits as N ^ oo of optimal flows across a typical A/ x M square 
define flows across M x M squares in the lattice I?. These are consistent 
as M increases and so define a flow on "from the infinite boundary to 
the infinite boundary." The constant 7(k, B) defined in Theorem 2 as a 
M — > CX3 limit could more abstractly be defined in terms of optimal costs 
in an appropriate problem involving such infinite-distance flows on 7?. See 
[3] for the details of such an approach in a different setting (traveling 
salesman problem through some specified proportion of random points). 

13. The abstract view above does not seem to help in the d-dimensional 
setting, but if instead the random network is "locally tree-like" (as with 
many models of random graph), it can be related to the cavity method [16] 
of statistical physics. That is, the n — > oo limit problem is an optimization 
problem on an infinite tree which can be tackled by setting up recursive 
equations exploiting the tree structure. Making such arguments rigorous is 
a major challenge; see [2] for nonrigorous methodology applied to optimal 
flow problems and leading to explicit numerical results. 

14. What kind of flow problems can be handled by the methodology of 
this paper, as outlined in Section 2.1? Intuitively, what seems important is 
the following: 

• We are studying an optimal global flow which minimizes a global cost 
function defined as a sum of local cost functions. 

• Constraints are local. 

• There is a stationary random environment determining costs and con- 
straints. 

In such a problem one can seek to modify Theorem 2 by replacing the specific 
definition of cm,b with its analog for different model assumptions. But we 
need the solution to have the property that flow volumes across different 
edges have the same order of magnitude, eliminating the "linear" case in 
comment 9 above. 
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